Nature Structural & Molecular Biology

Genes encoding novel protein sequences are a ubiquitous feature of genomes. They fuel molecular and cellular evolutionary innovations and frequently contribute to species-specific characteristics. We are now unravelling the processes by which they originate, including de novo from noncoding sequences and through extreme divergence, yet how much and what types of novel proteins evolve through each process is still unclear Does the mechanism of origination shape the structural and functional potential of the resulting proteins? Here, we conducted a broad computational investigation of genetic and protein novelty at the scale of the entire subphylum of Saccharomycotina yeasts. We detected more than 5,000 robust de novo genes across 332 species and compared them to more than 10,000 novel genes resulting from extreme sequence divergence, revealing two distinct modes of evolution of novelty. A remarkable 40% of de novo proteins are predicted to localize to mitochondria compared to only 15% of divergent, with the latter also being substantially longer and more disordered. A detailed analysis of conservatively predicted tertiary structures of novel proteins shows that "invention" of novel folds can happen through both processes but is more likely to occur de novo. We also illustrate cases of evolutionary "re-invention" of existing protein folds from non-coding sequences. Our work deepens our understanding of the origins and importance of novel proteins opening new directions for further structural and functional characterization.

Initiator-methionine excision is among the earliest protein modifications, yet its relationship to assembly geometry is unknown. Burial of the mature first residue was measured across 7,246 deposited human biological assemblies (22,291 chain-level observations; 1,191 proteins). Among 1,143 analyzable proteins, termini in MetAP-permissive penultimate-residue sequence classes were less often interface-engaged than termini in MetAP-nonpermissive classes (37.4% versus 47.4%; adjusted odds ratio 0.65, p = 7.2e-4). Curated processing annotations did not show a corresponding burial difference, and correlated residue properties preclude attributing the sequence-class association specifically to iMet removal. The analysis identified 264 interface-engaged MetAP-permissive candidates concentrated in cellular machines. In a fully recomputed conformer scan of deeply buried proteasome positions, modeled methionine accommodation was less favorable than at observed-methionine controls (median overlap -0.30 versus -1.12 angstrom, p = 0.0049), although most scoreable sites permitted a nonoverlapping placement. The census therefore reveals a graded structural constraint - not universal steric failure - and prioritizes complexes in which altered packing, assembly kinetics, lipidation or N-terminal methylation can be tested.

Life began with the emergence of a molecule that could replicate its own genetic material, a task plausibly mediated by an RNA-dependent RNA polymerase ribozyme. Here, we present the structure of such a polymerase ribozyme, bound to RNA substrates comprising the template, primer, and nucleoside triphosphate (NTP) analog. The structure reveals how directed evolution shaped flanking elements around a highly conserved catalytic core derived from the ancestral class I ligase ribozyme. Each element serves as a functional module, positioning the primer-template duplex and incoming NTP within the active site of the enzyme. This emergent domain organization is remarkably similar to the "right hand" configuration of polymerase proteins, suggesting a common functional form for copying nucleic acids, regardless of biopolymer catalyst.

Transcription factor Ets2 coalesce with the NF-{kappa}B pathway to regulate gene expression in specific signaling contexts. IKK2/{beta}-mediated phosphorylation events critically regulate the NF-{kappa}B pathway. However, any link between Ets2 and IKK2 remains elusive. Here we report Ets2 as a direct substrate of IKK2. In-vitro kinase assays using deletion constructs, high resolution MS-MS and site directed mutagenesis identified S295 as a prominent phosphorylation site distal to the DNA binding domain, substitution of which to phosphor-mimetic Glutamate triggers further phosphorylation of Ets2. MD simulations clearly indicate conformational constriction of the otherwise disordered N-terminal region and inhibition of DNA binding activity upon phosphorylation, which was further confirmed by Electrophoretic mobility shift assays. Our results uncover a phosphoregulatory connection between Ets2 and IKK2.

Human cytochrome P450 2C9 (CYP2C9) is a hepatic microsomal enzyme involved in the oxidative metabolism of clinically important drugs, but the structural organization of its oligomeric assemblies outside crystallographic packing environments remains poorly understood. Here, we report the cryo-EM structure of human CYP2C9 determined under aqueous, membrane-free conditions at 3.31 Angstrom resolution. The structure reveals a C2-symmetric hexameric assembly organized as a dimer of trimers. Individual protomers retain the conserved P450 fold and heme-binding architecture observed in previously reported crystal structures, indicating that assembly formation does not substantially perturb the catalytic core. The hexamer is stabilized by defined intra-trimer interfaces involving the N-terminal region and residues around Trp212 and Phe482, together with inter-trimer interfaces involving Leu71 and the 220-227 loop. These interfaces are distinct from the crystal packing contacts observed in CYP2C9 crystal structures, demonstrating that the assembly is not a simple recapitulation of crystallographic packing. Notably, the inter-trimer interface is located near the FG-loop-containing surface previously implicated in membrane association. This suggests that the observed hexamer may represent a membrane-free association of two trimers through membrane-related surfaces, whereas the trimeric arrangement itself may be compatible with membrane-associated organization. The structure therefore provides a framework for investigating how trimer formation, membrane interaction and local conformational changes in the FG-loop region may influence CYP2C9 function.

tRNAs are characterized by extensive chemical modifications that influence tRNA fate. N1-methyladenosine at position 58 (m1A58) is a widespread core tRNA modification linked to physiological and pathological processes. However, how m1A58 coordinate tRNA folding and processing to ensure translational efficiency in mammalian cells remains largely unknown. Using acute dTAG-mediated degradation and CRISPR-Cas9 knockout, we identified initiator methionine tRNA (tRNAiMet) as selectively vulnerable to m1A58 loss, lacking the isodecoder buffering observed for most other tRNA isoacceptors. NMR analysis of the tRNAiMet showed that m1A58 stabilizes D/T-loop interactions, consistent with a maturation-competent conformation. In vitro processing assays further demonstrated that m1A58 promotes RNase P-mediated 5'-leader removal and RNase Z-mediated 3'-trailer cleavage, while La/SSB protects accumulated precursors. Disrupting this checkpoint impaired the assembly of the eIF2-containing 43S pre-initiation complex and global protein synthesis, which was substantially rescued by adding m1A58-modified tRNAiMet. Acute TRMT6 degradation elicited temporally coordinated gene-expression responses involving proteostasis, transport and signaling. Together, these findings establish m1A58 as a conformational checkpoint coupling human initiator-tRNA maturation to translation initiation and stress responses.

Ammonia oxidation initiates nitrification and is closely linked to microbial N2O production. Ammonia monooxygenase (AMO) catalyzes the first and rate-limiting step of nitrification and is widespread across evolutionarily distinct ammonia-oxidizing archaea (AOA) and bacteria (AOB). The ocean is the largest biome for AOA and AOB, which have distinct ecological niches and markedly different sensitivities to nitrification inhibitors. However, the lack of archaeal AMO structures and inhibitor-bound AMO complexes has hindered mechanistic understanding of the architectural, catalytic, and inhibitory divergence between these two enzyme systems. Here, we report high-resolution cryo-electron microscopy (cryo-EM) structures of marine archaeal AMO captured in active and inactivated states within its native membrane environment, together with inhibitor-bound structures of estuarine bacterial AMO. Archaeal AMO forms an unexpected cup-shaped homotrimer composed of eight subunits per protomer and exhibits substantial architectural divergence from bacterial AMO. Integrated structural, biochemical, kinetic, and computational analyses reveal distinct periplasmic architectures, copper-center organization, and hydrophobic channels between archaeal and bacterial AMOs for ammonium acquisition, catalysis and inhibitor response. These findings provide a structural and mechanistic framework for understanding how archaeal and bacterial AMOs have diverged to distinct ammonia-oxidizing strategies and inhibitor susceptibilities across environmentally important ammonia oxidizers.

Canonical McrBC enzymes are nucleotide-powered, motor-driven endonucleases that bind and cleave modified bacteriophage DNA. Non-canonical McrBC homologs like LlaI and BsuMI are distinguished by a unique three-gene organization and the ability to target DNA site-specifically. Here, we report the atomic-resolution crystal structures of the DNA-binding module LlaI.R1 and AAA+ motor LlaI.R2 from the Lactococcus lactis LlaI anti-phage defense system. The crystallized LlaI.R2 hexamer traps two distinct active site conformations that correlate to different states of the nucleotide hydrolysis cycle and reveal that the organization of the critical catalytic machinery present in canonical McrB homologs is also conserved in non-canonical R2 proteins. Although canonical McrB homologs are strictly GTP-specific, we find that the R2 proteins from LlaI and BsuMI do not discriminate between different nucleotides, even when in complex with their respective R1 partners. Using mutagenesis, we define surfaces on the LlaI.R1 structure that are critical for DNA-binding and interaction with LlaI.R2. These observations support computational modelling of the assembled LlaI restriction system bound to DNA. Together, our data provide new insights into the evolution of nucleotide specificity in McrBC restriction complexes and the molecular mechanisms governing McrBC-catalyzed DNA translocation and cleavage.

mRNAs can form stable structures that need to be resolved to facilitate translation. During translation initiation in mammals, the scanning 48S complex requires the helicase activity of DHX29 to unwind stable mRNA structures that cannot be resolved by eIF4A. Here, we show that the yeast DHX29 homolog, Ylr419w (Dhx29), has a similar function during translation on elongating 80S ribosomes. Cryo-EM analyses show that the Dhx29 helicase module is positioned at the mRNA entry channel to engage mRNA, while its double-stranded RNA-binding domain (dsRBD) senses hairpin-forming mRNA in the ribosomal A-site. By selective ribosome profiling, we observed that Dhx29 is associated with transcripts that form RNA structures, such as stable tetraloops. Dhx29 mutants with perturbed helicase activity enrich 80S with hairpins in the A-site, as well as ribosome collisions, while a mutant lacking the N-terminal dsRBD sensor domain loses the specificity for such ribosomes. We thus propose that Dhx29 functions in translation elongation by resolving structured mRNA formed in the ribosomal A-site through its 3'-5' helicase activity and pulling on the mRNA from its 3' end.

Protein language models organize sequence and structure at scale, but a global representation of how proteins respond to mutation remains lacking. We present RegimeFormer, a large protein perturbation model coupled to RegimeAtlas, constructed by harmonizing and indexing 202,556,313 non-redundant protein sequences across the tree of life. A diversity-preserving one-million-protein subset provides the high-resolution training and inference layer, with 995,995 proteins yielding residue-level summaries across 407,048,356 residues and substitution-specific predictions available on demand. Across experimental deep mutational scanning, molecular benchmarks, structural confidence and evolutionary constraint, RegimeFormer identifies reproducible protein-level perturbation regimes that organize residue fragility, adaptability and predictive uncertainty. Regime conditioning improves substitution-specific prediction, with the largest relative gains under unseen-protein, unseen-family and low-homology evaluation. RegimeFormer-derived molecular priors further improve downstream transcriptomic and drug-response modelling. Together, RegimeFormer and RegimeAtlas provide a scalable framework for mapping, predicting and querying protein perturbation landscapes across global sequence space.

Gap junction (GJ) channels mediate direct intercellular communication by allowing the exchange of ions, metabolites, and signaling molecules between neighboring cells. Humans express 21 connexin (Cx) isoforms that can assemble into homotypic or heterotypic channels, creating a large potential interaction landscape that shapes tissue-specific communication networks. However, the rules governing which connexin isoforms can compatibly dock remain incompletely defined. Extracellular loop 2 (EL2) sequence features have been implicated in docking specificity and used to classify connexins into two canonical compatibility groups, K-N and H, but these assignments remain largely predictive. Most potential heterotypic connexin pairings have never been experimentally tested. This incomplete interaction map limits our ability to predict which connexin combinations can assemble, how isoform co-expression shapes intercellular communication, and how these relationships are altered or exploited in disease and engineered systems. Here, we used the FETCH (Flow Enabled Tracking of Connexosomes in HEK Cells) assay to evaluate docking compatibility across the complete human connexin family. To support family-wide compatibility mapping, we used literature-supported heterotypic interactions to define a data-driven FETCH score threshold for high-confidence interaction compatibility. Homotypic FETCH measurements varied substantially across the 21 connexin isoforms, with 15 producing mean scores above the empirical threshold. We then extended FETCH analysis to all 210 pairwise heterotypic isoform combinations. The resulting interaction landscape largely recapitulated expected motif-class relationships, including enrichment within the two canonical compatibility groups, but also identified neighboring-group interactions and unexpected cross-group pairings that represented clear exceptions to class-based predictions. Consistent with these findings, pairwise EL2 motif similarity was only modestly associated with threshold-based interaction classification, indicating that EL2 similarity alone was insufficient to predict compatibility outcomes. Together, these findings suggest that motif class provides a broad organizing framework for connexin compatibility, but that pairwise docking specificity also depends on yet-unresolved isoform-specific determinants that produce neighboring-group relationships and clear cross-group exceptions. Notably, Cx46, a lens Cx also associated with melanoma and breast cancers, emerged as a broadly permissive isoform capable of interacting with partners from both major compatibility groups and more than half of the connexin family. Together, these findings establish the first family-wide experimental atlas of human connexin docking compatibility, defining canonical interactions, previously unrecognized pairings, and exceptions to established compatibility rules. This atlas provides a foundation for defining the molecular determinants of connexin specificity, understanding how isoform diversity shapes intercellular communication, and designing gap junction channels with controlled docking behavior.

Apoptosis-inducing factor is a mitochondrial flavoprotein that links redox metabolism to mitochondrial homeostasis through its interaction with the disulfide relay protein CHCHD4. Although NADH-dependent AIF dimerization has been proposed as the activated state mediating CHCHD4 engagement, whether it is strictly required for productive AIF-CHCHD4 function remains unclear. Here, combining cellular, biochemical and biophysical approaches, we show that disruption of the AIF dimer interface compromises oxidative phosphorylation, respiratory-chain organization and CHCHD4-dependent mitochondrial homeostasis, yet preserves partial AIF function. Our data reveal that the AIF-CHCHD4 system operates as a conformational dynamic redox module in which distinct AIF oligomeric and redox states sustain CHCHD4 activity with different efficiencies. Mechanistically, dimerization is coupled to NADH-dependent conformational changes that regulate coenzyme binding, charge-transfer complex stabilization and catalytic efficiency. In turn, CHCHD4 binding remodels AIF conformational and redox properties, partially compensating for defects in dimer stabilization or redox coupling. Consistently, a peptide derived from the CHCHD4 N-terminus partially restores redox function in a pathogenic AIF variant defective in dimer stabilization, supporting partner-assisted allosteric regulation as a potential therapeutic strategy.

UCH-L1 is a monomeric deubiquitinating enzyme whose native structure embeds a shallow $5_2$ knot located near the N-terminus, placing the knotted topology in direct proximity to both the substrate-binding pocket and the catalytic site. While our previous work established that N-terminal integrity is critical for catalytic activity, the energetic cost of unknotting and its structural consequences remained unquantified. Here, we combine steered molecular dynamics with an umbrella sampling scheme to generate topologically modified variants of UCH-L1 and, for the first time, reconstruct the free-energy profile of UCH-L1 unknotting. The potential of mean force reveals a steep energetic barrier to knot disruption, consistent with knotting being a late, rate-limiting folding step that is effectively locked in once the native structure is established. Long unbiased MD simulations of fully unknotted variants in both apo and holo states show that knot removal increases local flexibility at the N-terminus without inducing significant global structural destabilization. Binding energy calculations indicate that the unknotted variant binds to ubiquitin less tightly than the wild-type ($\sim$-62~vs~$\sim$-76~kcal/mol), suggesting that topological integrity contributes to substrate affinity. Together, these results show that the $5_2$ knot in UCH-L1 is not a passive structural feature but a functional element that fine-tunes folding kinetics and contributes to substrate binding efficiency.

The 5S ribosomal RNA gene is present in the human genome not once but in ~80 copies, arranged head to tail in a single array of ribosomal DNA on chromosome 1 -one of the most repetitive and least explored regions of the genome. Its product is one of the four RNAs in every ribosome and, when ribosome assembly fails, it activates the tumour suppressor p53. Whether these copies vary in sequence between people, and whether such variation has physiological or pathological consequences, is unknown. Using telomere-to-telomere genome assemblies, whole-genome sequences from ~490 000 UK Biobank participants, and ~940 GTEx transcriptomes, we find that every person carries copies bearing substitutions or indels, and that ~10% of people express such variant 5S rRNA. Mutating every position of the gene in vitro, we find that variants blocking incorporation into the ribosome map to the uL5/uL18 interface and activate p53. Remarkably, these same variants are depleted from human populations: selection has acted on the step that p53 monitors. Ribosomal DNA is thus a functional source of human genetic variation, long invisible to genome-wide analysis and shaped by the p53 pathway it controls.

Linker chemistry and conformation are central determinants of PROTAC activity, shaping ternary-complex geometry, cooperativity, target-lysine presentation and cellular permeability. Existing linker generators often lack explicit control over linker flexibility, require predefined attachment sites and linker lengths, or produce structures that demand substantial geometric correction, limiting their utility in practical PROTAC design. Here we introduce FlexiTAC, a Bayesian flow network that jointly generates linker atom types and coordinates from the warhead and E3-ligase-ligand contexts. We also assemble PROTAC-3D, a quality-controlled collection of 63,554 component-resolved PROTAC structures for model training, and PROTAC-Bench, which covers molecular quality, fragment preservation, geometric fidelity, conformational stability, fragment awareness, rediscovery and sampling efficiency. Compared to the best 3D baseline models, FlexiTAC improves validity by 12.0-12.7% and achieves the highest PoseBusters pass rate of 79.5%-80.0%. A differentiable guidance module shifted generated linkers along a conformational ensemble-derived rigidity axis without retraining the generator. In silico case studies further show that the model can accept crystal-derived, redocked or predicted structural inputs. Together, FlexiTAC, PROTAC-3D and PROTAC-Bench establish an integrated and reproducible framework for data-driven PROTAC linker design, combining controllable structure-conditioned generation with standardized training data and evaluation protocols. This framework expands the linker chemical and conformational space accessible to computational exploration, provides a foundation for future method development and enables the systematic generation of structure-conditioned linker designs with tunable conformational flexibility.

O-linked {beta}-N-acetylglucosamine (O-GlcNAc) transferase (OGT) is the sole writer for intracellular O-GlcNAcylation. It catalyzes O-GlcNAcylation of thousands of protein substrates, but relatively less is known about the post-translational modifications that occur on OGT itself. Herein, we demonstrate that OGT is S-palmitoylated at Cys-472 and Cys-477, which is mediated by the S-acyltransferase Zinc Finger DHHC-Type Palmitoyl transferase 14 (zDHHC14) and removed by acyl protein thioesterase 2 (APT2). S-Palmitoylation stabilizes OGT by shunting it away from the lysosomal chaperone-mediated autophagy (CMA) pathway, as S-palmitoylation decreases the interaction between OGT and heat shock cognate 70 kDa protein (HSC70), the CMA chaperone. Via label-free quantitative mass spectrometry, we find that S- palmitoylation elevates the affinity between OGT and protein phosphatase 1 catalytic subunit gamma (PPP1CC), but not PPP1CB. We further demonstrate that S-palmitoylation of OGT augments binding with Yes-associated protein-1 (YAP), a protein that associates with PPP1CC, and subsequently enhances YAP O-GlcNAcylation. Our work unearths S-palmitoylation of OGT and CMA-mediated degradation of lysosomal OGT, the orchestration of which finetunes the activity of key OGT complexes, such as OGT-PPP1CC, and contributes to OGT substrate selectivity.

Oil palm (Elaeis) is the primary source of global vegetable oil. Interspecific hybrids of Elaeis exhibit pronounced heterosis by integrating two distinct subgenomes into a single nucleus, effectively combining the high yield of African oil palm (E. guineensis) with the high unsaturated fatty acid content and disease resistance of American oil palm (E. oleifera). However, the genetic basis underlying heterosis is still unclear. Here, we combine phased genome assembly, comparative genomics, evolutionary genomics and haplotype-aware transcriptomics to unravel the genetic architecture of heterosis of hybrid oil palm. We assemble the highly heterozygous F1 genome ('Reyou 40', 3.75% heterozygosity) into a complete 1.73 Gb T2T haplotype (HapG) and a 1.84 Gb near-T2T haplotype (HapO with17 gaps). Despite 91.56% sequence identity, HapG and HapO diverged in LTR-RT occurrence and PAV affected genes, showing complementary biases in lipid metabolism and stress responses, respectively. Evolutionary genomics revealed that ancient WGDs preserved the palm family. Whereas lineage-specific lipid-related gene expansions in oil palm. Six ancient introgressed regions (~64 Mb) in HapG were reshaped by transposable elements and tandem duplication, showing an enrichment of genes related to resistance and lipid metabolism. Transcriptomically, 82.2% of allelic gene pairs maintained balanced expression, accompanied by parental functional complementarity and dosage buffering, revealing a potential regulatory basis for coordinating parental genetic differences in the hybrid genome. These haplotype-resolved genomic resources offer vital targets for understanding heterosis and accelerating oil palm molecular breeding.

IgM is the most ancient antibody isotype, playing an important role in both circulatory and mucosal immune responses across vertebrates, yet structural characterization of its polymeric forms is limited outside of mammals. Here, we report the cryo-electron microscopy structure of mallard duck secretory (S) IgM at 3.37-[A] resolution. The structure revealed a pentameric core globally similar to human SIgM, supporting the view that pentameric IgM is subject to strong evolutionary constraints. However, compared to mammalian structures, we observed species-specific differences at molecular interfaces. Surface plasmon resonance binding assays characterizing secretory component (SC)-IgM interactions supported structural observations and, when compared to IgA binding, revealed isotype-specific contributions from the avian SC N-terminal extension. Together, these findings establish a comparative structural framework for polymeric IgM across vertebrates and provide insight into how avian SIgM-specific features may support mucosal immunity in birds.

Disulfide-rich peptides possess exceptional potency and selectivity but are often limited by the instability and synthetic challenges associated with native disulfide bonds. Here, we report the design, synthesis, pharmacological evaluation, and structural characterisation of triazole-based peptidomimetics of the -GI conotoxin, a selective antagonist of the muscle-type nicotinic acetylcholine receptor (nAChR). A series of 1,4- and 1,5-disubstituted triazole analogues were prepared entirely on resin using CuAAC and RuAAC chemistry to replace the native Cys3/13 disulfide bridge. Functional evaluation against human muscle nAChRs revealed that 1,5-triazole analogues retained low-nanomolar potency, with the lead mimetic exhibiting activity comparable to native -GI. Cryo-electron microscopy of the lead compound bound to the muscle-type nAChR provided the first structure of a disulfide-isostere peptidomimetic in complex with a membrane receptor. The structure demonstrates that the 1,5-triazole reproduces the native peptide fold with high fidelity while contributing receptor-facing interactions not available to the native disulfide bridge. Molecular dynamics simulations further revealed conserved hydration networks and similar conformational sampling between the native peptide and lead mimetic. Together, these findings establish triazoles as effective disulfide surrogates and provide a structural framework for the rational design of stabilised conotoxin therapeutics.

A developmental increase in cerebrospinal fluid (CSF) production during development is essential for neuronal growth and ventricular expansion. A key regulator of CSF production is the specialized sensory multicilia of the choroid plexus (ChP), which mediate non-canonical Sonic hedgehog (Shh) signaling to suppress water channel and ion transporter expression, thereby limiting CSF production. ChP multicilia progressively shortens during development, attenuating Shh signaling and promoting CSF production. Here, we identify miR-34/449 miRNAs as essential regulators of ChP multiciliogenesis. Whereas mutations in canonical ciliogenesis genes elevate CSF production and contribute to hydrocephaly, deletion of miR-34/449 reduces CSF volume and causes microcephaly. Loss of miR-34/449 miRNAs causes excessive basal body amplification, defective basal body docking, and failure of developmental multiciliary shortening. Consequently, miR-34/449-deficient ChP cilia remain abnormally long and fail to attenuate Shh signaling, resulting in sustained repression of water channel and ion transporter expression and reduced CSF production. Mechanistically, miR-34/449 miRNAs directly target Gmnc, a master transcriptional regulator of multiciliogenesis, to restrain basal body amplification and promote basal body docking. Together, our findings identify miR-34/449 miRNAs as critical regulators of ChP multiciliogenesis and establish the developmental remodeling of ChP multicilia as a mechanism to couple Shh signaling dynamics to developmental control of CSF production.